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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

TECHNICAL MEMORANDUM 1288 

THE DIFFUSION OF A HOT AIR JET IN AIR IN MOTION* 
By W. Szablewski 

PART II. THE FLOW FIELD IN THE TRANSITION ZONE 

SUMMARY 


The turbulent diffusion of a hot air jet in air can be divided into 
two zones, the core and the transition zone. The first part of this 
study (reference l) deals with the flow field in the core, the second 
(the present report) with that in the transitional zone. 

Part A of the present report is limited to small temperature differ- 
ences. The decrease in the velocity and the temperature along the jet 
axis, the breadth of the mixing region, as well as the asymptotic distri- 
bution functions, are determined. 

The empirical constant K, a measure for the mixing length, appearing 
in the theory, follows closely a value of 0.010 for asymptotic conditions. 

Experimental data ,are available only for the case of outside air at 
rest. The comparison with theory indicates that the asymptotic distri- 
butions are satisfactorily reproduced with exception of the boundary zone. 
Velocity and temperature drop are very closely reproduced by the computed 
functions up to a point near the boundary of the core, while a last short 
fraction of the region adjacent to the boundary of the core is not covered 
by the theory. 

Part B of the present report deals with greater temperature differences. 
The breadth of the mixing region, as well as the velocity and temperature 
drop along the jet axis, is calculated. The theory is then compared with 
Pabst*s measurements (reference 2). The ratio of interchange of temperature 
and velocity yielded a factor E = 2. Considering the friction loss at the 
nozzle wall, the agreement between the theoretical and the experimental 
velocity decrease along the jet axis can be regarded as satisfactory. The 
temperature measurement along the jet axis appears to be faulty. 


"Die Ausbreitung eines Heissluftstrahles in Bewegter Luft." 
GDC/2if60, September 19^-6. 
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A: SMALL TEMPERATURE DIFFERENCES BETWEEN JET 

AND SURROUNDING MEDIUM 
I. Method and Results 

1. In the first part of the investigation (reference l) the flow 
field in the core was computed. 



The investigation included the variation of the curves hounding the 
mixing zones of velocity and temperature as well as the aspect of the 
velocity and temperature distribution functions over the mixing zones. 

In the second part, (the present report), the flow field in the 
zone adjoining the core is investigated. 



This zone, which in the asymptote in the so-called axially symmetrical 
jet diffusion is characterized by the affinity of the flow processes in 
the cross sections of the jet, is termed the transition zone. 
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Variation of mixing width of the velocity and the temperature over the 
nozzle spacing 


The following relations are involved: 
bp the breadth of the mixing region of the velocity 

bp the breadth of the mixing region of the temperature 

In addition 

•a o temperature rise of the discharging jet 

■d temperature rise of the jet on jet axis 
u Q velocity of the discharging jet 

Up velocity of the surrounding medium 

u^ velocity of the jet on jet axis 



Variation of jet velocity and jet temperature along the jet axis 
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2. In the extension of the theory of free turbulence to gases of 
widely variable densities, the following turbulent exchange quantities 
had been obtained in part I (reference l). 


M = El 2 


du dp 
dy dy 


for turbulent diffusion (la) 


T = V 


I! (f + Et f 


Q = El 2 


du 

d(pcpT) 

3y 

Sy 


for turbulent shearing stress (lb) 
for turbulent heat conduction (lc) 


On the assumption of constant pressure for the diffusion of a hot 
air jet in air in motion (axially symmetrical case) the motion equations 
were then obtained: 

equation of continuity of mass 


d(rpu) + d(rpv) = Ee ( x ) d_ /rdp\ 
dx dr dr \dr ) 

equation of continuity of momentum 

d(rpuu) , d(rp-\ni) , / -d u ~-8p 

- - - — - + — = e(x)2L_ rp — + Eru— 

dx dr dr \ dr dr 

equation of continuity of the heat (energy principle) 


d(rpuT) d(rpvT) 

+ 

dx dr 


Ee (x)|- 
or 


r 8(pT) 

dr 


with the apparent kinematic viscosity 


€ (x) = Kbp 


Umax ~ u min 


Owing to the continuity of mass, it further yields: 
momentum 


(2a) 


(2b) 


(2c) 


drpu^u - ujj drpv^u - up'j ^ ^ 


dx 


dr 


dr 



II I III I 
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Heat 

d(p|) 

8r 

(•8 = temperature rise of jet.) 

Integration with respect to r gives then for the transition zone: 
momentum 

(3a) 


rpv(u - u£) - 



ax. c U) r fp%^) + K( 5 - Ul )| 


+ = Ee( x )Z- 

dx dr dr 


heat 


rpv-g - r~ / pu-3r dr = Ee(x)r^— (pd) 
ox J r dr 

with e(x) = Kb-^^u^ - u-j^ < 

Momentum and heat law can he physically interpreted as follows: 


(3b) 



On marking off a control area in the mixing zone in the manner 
indicated above, the momentum theorem states that the loss of momentum, 
which the flow experiences on its passage through the control area, finds 
its equivalent in the turbulent shearing stress 
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at point r. Correspondingly the heat balance goes through the control 
surface with consideration of the heat convection only at point r when 

the turbulent heat conduction Q = Es(x)^— (p 6 ) is included. 

dr 

For r = 0 a further integration is possible which gives the 
formulas of the conservation of momentum 




dr = 


P 0 U 0 (uo - 



and of heat> 



pudr dr = PqUqtS 0 



( 1*0 




With the equation of state for perfect gases at constant pressure as 
basis 


p T = const ( 5 ) 

p can be replaced by in the above equations 

- _ const _ const 1 
P ” T Tq 1 +($/T{) 


where is the temperature of the surrounding medium. 


3. For the investigation of the transition zone by means of the 
momentum and heat equations as well as of the equations of conservation., 
the premise 


u - ui u A (x) - Up 

= — 9(r,x) 

u o ' U 1 u o " U 1 


( 6 a) 


■dho = ^(r,x) ( 6 b) 

^o 

is made, where u A (x) is the jet velocity on the jet axis and t 3 a (x) 
is the temperature rise on the jet axis. 
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This formula, carried in the equations, gives (if a practicable 
assumption can be made for the variation of cp(r,x) and ty(r,x)) four 
equations for the four unknowns: 

? A (x), $ A (x), b-^x), and b 2 (x) 


The transverse component v still appearing in the equations is 
determined from the equation of continuity of the mass 

v = ^ r dr (7) 

r Jo 


The next problem is to make a practicable premise for cp ' and 
These functions occur in the integrands of the integrals appearing in 
the equations and in the expressions defining the shearing stress and 
heat conduction at point r, whereby a mean value may be chosen for r. 

In connection with the similarity of flow in the jet cross sections 
resulting for the asymptote it is then logical to put 


cp = cp(T]*) with 

II 

* 

p- 

(8a) 

\|r = t|/(t}) with 

r 

^ =bi 

(8b) 


This theorem actually seems to be confirmed by experimental results which 
indicate far reaching similarity of flow profiles in the transition zone 
(reference 2). 

Recommended for the profile form is the' distribution function resulting 
theoretically for the asymptote: for the case of air in motion (see 

section III) It is the function e"*^ OT ^ , where the parameter a, defining 
the width of the distribution, contains characteristics of the asymptotic 
flow. The actual calculation by this method presents some difficulties, 
and besides the distribution function obtained by the differential equation 
does not approximate the experimental distributions in the desired measure. 
(See fig. 9.) 

Incidentally it seems timely to comment upon a manifest inadequacy 
of the theory, which may be attributable to a not entirely equivalent 
formulation in Prandtl's concept of the mixing length. According to 
Prandtl’s hypothesis of the mechanism of turbulent mixing, the 
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turbulent exchange is effectuated by turbulence balls, which for the 
duration of their own existence transfer their properties invariantly 
from one layer to the other. Accordingly the momentary fluctuation u* 
with the momentary mixing length 2 should be put as 

u* = u(y + 2) - u(y) = + . . . 

In the original mathematical wording u' was put equal to 2 Si 

dy 

which, in general, represents a permissible approximation considering 
the smallness of 2, but at particular points such as in the profile 

center, for example, where = 0 and therefore u* should be 

dy 

put u' = 0, the results are inadequate, although fluctuations are 
certain to occur. In consequence, Prandtl (reference 3) suggested the 
average value formation 

e ^ u max “ u min| 

(wherein K = empirical constant, b = breadth of mixing region of the 
velocity) for the apparent kinematic viscosity 

7 2 du du 

£ — L 

dy dy 

in the original draft. This theorem has then no longer the differential 
character of the original formulation and avoids the said inadequacy; 
however, it then results in inadequacies at the edge of the profile where 
the fluctuations and, with them, the apparent kinematic viscosity cancel 
out, whereas by the new theorem an amount of e constant over the entire 
width of the profile is involved. On the other hand, the new theorem has 
the implicit advantage of being substantially simpler mathematically and 
hence of simplifying the calculations considerably. 

A very practical representation of the distributions, which gives 
a very good approximation, is presented by the function 

^1 - T]* 3 / 2 ) 2 and ^1 - t] 3//2 j 2 

already used by other workers (reference k). 
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Thus 

y!f = (l - ti 3/ 2 ^ 2 (9a) 

q> = (i - n* 3//2 ) = |jl - (Ti-h 2 /bi) 3 / 2 J 2 (9b) 

with t] = r^b 2 . To simplify the calculation for b 2 ^b^ for the case of 

outside air in motion fu^ ^ 0 ), the ratio of the widths resulting from 
the theory of the asymptotic distribution functions for outside air in 
motion is approximately put as 


b 2 /k>i = y/E" (10) 

(See section III.) 

The case of quiescent outside air ( U 1 = o) is treated later. 

k. Part A of this report is restricted to small temperature differ- 
ences. This case is amenable to calculation and represents the conditions 
accompanying any temperature differences in first approximation. 

The velocity field is computed by the momentum equations. 

The equation of conservation of the momentum gives the breadth of 
the mixing region of the velocity b^ as a function of the axial 
velocity u^ 



( 11 ) 


with the constants dj = 0.133 and e 1 = 0 . 257 * and the subsequent 
calculation gives the axial velocity u^ as the describing variable 
of the flow field rather than the nozzle distance x. The momentum 
equation is specialized to r = r D . 

The use of the foregoing result gives a linear differential equation 
of the first order for the nozzle distance x as function of the velocity 
on the jet axis u^. 
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When the integration constant is fixed by the postulate 
for the boundary of the core x^. 


U A ~ U 1 
u o - U 1 


= 1 


with 



( 12 ) 


and the constants 
a = 0.0275, 0 = 0.05^9, 


7 = 0.0375, 


H = 1.9259 


The length of core x^/r D is to be taken from the theory of the core; 
K is an empirical constant. 


In the asymptote the formulas are obtained 

*1 Iro X&f 2 *! i 7- 

1 WJ K - ui f/ 2 

[^o - u l) 

with the constant A^ = 1.972 
and 


(13) 


2L * l ( *0 \ h\ l/2 A 1 

r Q K\uo - up) \u 0/ l 2 


u^Up\ 3/ 2 

u o - Up) 


(14) 


with A2 = 0.1002 
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With the nozzle distance as independent variable the asymptotic 
formulas read 

/ Up - up j 1 / 3 

b l/ r o ~ K 1//3 - ---7T - BlfcAoV^ 3 (15) 




with B]_ = 4.246 and 


U A - U 1 ^ 1 

Uo ~ U1 k 2/3 


with B 2 = 0 . 216 . 


(ui/uq ) 1 / 3 

u o - u l\ 2/3 
U o ) 


B2 


( X / r o) 


273 


( 16 ) 


In figures 1, 2, and 6 the functions of the velocity mixing field 
are shown for the parameter values 


Uq - ui 

- - - = 1.0, 0.75, and 0.5 

u o 

The temperature field is computed by means of the heat equations. 

By the equation of conservation of heat, the breadth of the mixing 
region of the temperature bg is 


b 2/ r o = 



ui 


d 2 + q iL 


1/2 


with dg = 0.0786 and e^ = 0 . 257 . 


(17) 


The heat formula contains the empirical constant E, the ratio of 
interchange of temperature and velocity; according to the measurements 
by Pabst (reference 2) and others it is equal to 2. 


From the heat formula a Bernoulli differential equation is obtained 

for d fa as a function of Unfortunately the quadrature cannot 

A/ o u o “ u l 

be carried out. 
tion -9 . /9 ~ = 1 


V>o - 


Fixing the integration constant by the initial condi - 

■ u* - u n , , . 

A = 1, the result is 


for 


U - Ui 

o -L 
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with the functions 



and the constants 

C = -0.6111 D = 0.4848 

M.-L = 0.0241 H 2 = 0.0522 n 3 = O.O876 = 0.2571 

v ± = 0.0066 v 2 = 0.0253 = 0.0196 

Vl, = °- 1 335 v 5 - ^ v 6 ^ v 7 - >V 

In the asymptote the formulas 


h 2 / r o ~ \/2 *ifr 0 


(19) 
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are obtained in agreement with the mathematical assumption 

^ 2 / b i = vE-V? 


In addition 


® A^o ~ 



( 20 ) 


In figures 3> 5> and 6 the functions of the temperature mixing 

field are represented for the parameter values 


uu - u- 


1 _ 


uo 


1.0 


0.5 


The integral occurring in formula (l8) was obtained by numerical inte- 
gration. 

The case of quiescent outside air fu^ = 0) presents a singular 

behavior as evidenced by the fact that the breadth of the mixing region 
in the asymptote is represented by a linear function bj(x) ~ x, while 
with outside air in motion bj_(x) ~ x^/3, it is found further that the 

ratio of the asymptotic mixing width obtained for outside air in motion 
is not applicable here. The theory of the asymptotic distribution 
functions produces, in this instance, an impracticable result (t> 2 1 ^ 1 — * 00 ) 
The asymptotic ratio 


b 2/ b l = 1# 33 for u-j_ = 0 ( 2 1) 

is obtained from the law of conservation of heat by an approximation 
method. 


The functions of the turbulent diffusion for the singular case 
Uj = 0 are as follows: 


velocity field 



( 22 ) 
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with F x = 2. 737 



with F 2 = O.13V7. 

With the nozzle distance xyt Q as independent variable 




with G 1 = 2.737, G 2 = 20.321 




( 23 ) 


( 24 ) 


( 25 ) 


( 26 ) 



( 27 ) 


with 


% = 3.245 


h 2 = 1.123 
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For infinitely great nozzle distances, the formulas are 


b l/ r o ~ KG^ro ( 28 ) 


and 


with Kj_ = 3.644 



A> 



(29) 


(30) 


(31) 


with K 2 = 0.793. 

All formulas still exhibit the empirical constant K. It is 
embodied in the apparent kinematic viscosity 


e(x) = Kb-j^UA - u-jJ 


and represents a measure for z/b, the ratio of mixing length to mixing 
width, which is regarded as constant for the individual jet sections; 

K and Z/b are to be considered functions of the characteristic 
length x/r Q . But in view of the far-reaching similarity of the profiles, 

the dependence on x/r Q throughout almost the entire transition zone is 
expected to be slight, as the measurements seem to confirm. For example, 
Tollmien's investigations (reference 5) at the plane jet boundary (these 
conditions prevail in the immediate vicinity of the nozzle mouth.' ) 
give Z/b = 0.068; for the axially symmetrical jet diffusion (the 
conditions encountered at very great distance from the mouth of the 
nozzle), he obtained Z/b = 0.073 (where, for reasons of continuity 
at transition from the core to 'the transition zone, b in the latter, 
as done here, is to be put equal to the jet radius). 


Also of interest are the distribution functions (section III). 


already pointed out, the empirically obtained function 


(l - t,3/2) 2 


As 

gives 
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a very close approximation of . the experimental distributions. The theory 
gives the following asymptotic distributions for u^ f 0: 

Velocity distribution 


<P 


= e'^o 1 !) 


2 



(32) 


whereby with 


U A - ui 

Uq - Ui 



b-L z b M xl/3 


CT 


O 


-i JTk 

V 6 K u. 


^1 


As reflected by the previously obtained results concerning the structure 
of the mixing zones 


a Q = 1.811 


(33) 


for q = rj b^. 

Utilizing in correspondence with the asymptotic behavior of the 
mixing width ( b-j^ ~ b ro X x 1 /^) the coordinate 



with 


cp = e-^o ^*) 2 


CT * 

o 


0.427 


#73 





(34) 


In the case of quiescent outside air 
velocity distribution function is 


(“i - °) 


the asymptotic 



(35) 


■■■in mi h lima i mu hi an mi mi 


ill II ■ I l III 
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with. 



and 

ct d = (i/k)x 0.0784 

For the temperature distribution the important relation 

* = cp 1 /® and ^ = qp 1 / 2 (36) 


is obtained. This result is already indicated in Reichardt's report 
(reference 6). 

The asymptotic distribution functions are represented in figures 7 
and 8. 


Now, the theory is compared with the experimental results. The 
former contains two empirical constants E and K, E being equal to 2. 

In the first part of the investigation (reference l) the constant K 
for the core had been determined on the basis of Tollmien' s study at the 
plane jet boundary (reference 5). The result was K = 0.0106 on the basis 
of the heat equation, but the thus defined value is somewhat uncertain 
owing to the doubtful flow losses as a result of the friction at the 
nozzle wall and the assumption — valid exact only for small velocity and 
temperature difference — that the breadth of the mixing regions of 
temperature and velocity which served as basis of the calculations 

act as : 1. 

Quantity K is determined next for asymptotic conditions. Tollmien 
(reference 5) obtained b . = 0.2l4x for the diffusion of a jet issuing 
from point source, as against bp = K x 20.321 X x according to the 
calculations by equation (28). The comparison gives K = 0.0105. When 
the determination of K is based upon experimental results of other 
structures of the diffusion field, such as the gradient of the asymptotic 
distribution function or the decrease of velocity and temperature along 
the jet axis, the foregoing value of K is almost exactly reproduced. 

The calculations have been based on the value 


K = 0.010 


(37) 
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For asymptotic conditions, no incidental flow loss due to friction 
at the nozzle wall needs to be considered. (The flow loss can be 
allowed for by introducing an effective nozzle radius (see part I, 
(reference l) ) ; for infinitely great distance from the nozzle, however, 
the size of the nozzle radius is of no influence.) The ratio of the 
breadth of the mixing region of temperature and velocity serving as a 
basis of the calculations is rigorously valid for the asymptote. The 
above value of K should therefore represent a safe value for asymptotic 
conditions. 

Experimental data for small temperature differences as treated here 
are few and limited to the case of outside air at rest ^u^ = 0 j. For 

greater temperature differences (to be discussed in part B of this 
report), Pabst’s comprehensive measurements (reference 2) are available. 
For the asymptotic velocity distribution Reichardt’s measurements 
(reference j) are available. The comparison (fig. 9) shows practicable 
agreement up to the boundary zone where the divergence is fairly great. 
This difference is, as stated above, attributable to the nature of the 
calculation method. 

The decrease of velocity and temperature along the jet axis was 
measured by Ruden (reference 8). Unfortunately the published report 
dealt only with the test curves without giving the test points or any 
further details of the measurements. The comparison reproduced in 
figures 10 and 11 indicates very good agreement between the test curves 
and the computed functions up to a short transition arc. 

This transition arc still defies interpretation at the present state 
of turbulence research. The calculation by integral formulas contains as 
essential premise the assumption of profile similarity for the individual 
jet sections. The result is a hyperbolic variation of the distance 
function up to the boundary of the core, while the experiment and also 
the differential equations (2) of the mixing process indicate zero tangent 
in the boundary of the core. It is in this short transition arc that the 
transition from the asymptotic profile form to that of the plane jet 
boundary is largely effectuated. The profiles throughout the entire core 
are very similar to the profile forming in direct , proximity of the nozzle 
mouth, which corresponds to the mixing of two plane jets (compare part I 
(reference l)). For the calculation of the transition arc, a return to 
the differential equation (2) would be necessary, while the determination 
of the breadth of the mixing region involved in the apparent kinematic 
viscosity e(x) = Kb-^u^ - u-j_j, would call for the conservation formula. 

Even so, the variation of K within the arc would remain, whereby it is 
to be noted that the assumption of similar profiles is equivalent to 
assuming a constant K, so that the existing dependence of K on x/r Q 
would essentially have to occur inside the arc. But every hypothesis 
for this dependence of K and of l / b (assumed as constant in the 
individual jet sections) fails at present. 
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According to the experiment, the transition arc joining on the 
common core boundary terminates in the hyperbolic branch of the distance 
function; the velocity decrease in the arc is at first very slow compared 
to the temperature. In the comparison (figs. 10 and ll) the starting - 
point of the computed hyperbolic branch was chosen accordingly. 

In the calculation, the integration constant was so defined that 
— u. - u-, 

0 = 1 for — — = 1, i.e., that the hyperbolic branches of the 

' U o ~ U 1 

velocity and temperature decrease started at a common boundary. In these . 
conditions, a displacement of the experimental correlation is to be 
expected relative to the theoretical correlation of the related velocity 
and temperature, which would have to disappear for greater nozzle dis- 
tances. This displacement is manifest -in the comparison, figure 12. 


II. Calculation of the Variations of the Mixing Width 
and of the Axial Functions 


Momentum equation (3a) 

rpv^u - u-j^ - Z^ 1 pu(u - up)r dr = s(x)r 


* »p - 


Heat equation (3b) 


where 


-v 1 -j 2 _ 

rpv8 - — / pu^r dr = Ee(x)r — (p^) 

dx , L dr 


e(x) = Kbp (u A (x) - u-jJ 


Equations of conservation (4a) and (4b) 


pu^i - u : Jr dr = PqUo^o - up)r 0 2 y/; 


n b 2 , 

puflr dr = Pq^o r o 2 / 2 
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In correspondence with the equation p X T = constant p is replaced by 


const 1 
T x 1 + fr/Tp 


The transverse component is obtained from the equation of continuity 


d(ru) + d(rv) = Q 
dx dr 


1 / V du . 
v = - — / — r dr 

r Jo 8x 


Hence, for the momentum: 

du 


dx 


r dr 


(u - Up\ 

V J 1 + tf/Tp 


= Kb- 


. ( U A " u l) : 


1 + ti / T n dr 



r dr 


( 38 a) 


and 


J 1 a /r. 


( 5 -..yi. r ar = V (^,- , ^ r 0 2 /2 


/o 1 + V T 1 1 + ^o/ T l 


( 38 b) 


for the heat: 


du . \ 'A 
-r dr 


dx } 1 + til'll dx 


J 2 _ ti 

u = r dr 


EKb 1 (uA - up) r 2_ 


1 + tij’. 


dr 1 + ti j Tp 


(39a) 
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and 


- i 

/ - u z=~. — r 

l/0 1 + V T 1 


•a 


to = Uq 


1 + 




(39b) 


If we limit ourselves to small temperature differences 
we get: momentum 


(V T i)« 1 > 





and 



(1+Oa) 


( 4ot>) 


heat 



and 


’ 2 / 


\ 


-/N h r to 


Posing the formulas (6) and (8) 

■t (ti) 


ik 

r» - — r 



(4la) 


(4lb) 


with 


n = r/K 2 (x) 
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the transformation formulas read: 


n = i/bg, r = T)b 2 ( J+2) 


dr = bg dr} 

( pMt - (V®),, - (°/ ;h i) x i 


(a/a^) x = (a/a»i) x 


Corre spondingly 

n - U 1 = - U 1 cp(|) 

u o “ U 1 u o " U 1 

with 

6 = r/ b l = t) bg^bj 


giving, after effecting the transformation 
momentum 



(Accents denote derivatives with respect to x) 




2b 
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Here (9) must "be Inserted: 


t(*l) = (l - il 3 / 2 ) 2 . 


^ +* . e - [? b ^ 3/2 ) ° ;v bl/b 




bi ho = = i 


(a) Velocity Field 

I. Equation of conservation of momentum (43t>) 


- Ul\ 

Uq - Ul/u 

A “ U 1 


u o V u 

o 

i 


= 1 


The evaluation of the integral gives 

J 5 q» 2 e dl = f 1 (l - |3/2) 4 | as . | |2 - 8 6 


7/2.5 o 13/2 8 

+ - I - Al + 1 6 
5 13 8 


i ^ 4 7/2 i ^ 

(Pi d| . (1 - I 3 / 2 ) I di * | I - = I + i I 


hence 


cp 2 ! d| = 


243 

3^0, 


/. 


0 « 4 ‘-i5 
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with dj = O.I 3352 and e 1 = 0.2571^. 


In addition 
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Specialize to 5 = l/2. 
Evaluation of integral 
( 1 ) 


1/2 . 
<PS 


10 


21 

160 


~^2 = 0.0807k 
28 


( 2 ) 


9 


5=1/2 ( 


> A - f 3 / 2 \ 


j |=l/2 8 2 


- I - |iP = 0.4179 


^ 46 dlj <Pjy 2 = 0.03374 


' 1/2 *£ -2 


J g i 2 al - / 3 (l 2 - 5 1 / 2 ) i 2 ai , |! 5 - || 7/2 

|2 dE ' lio - |g 1/2 =-°.°57 01 

^/ 1/2 S ^ dE j V 2 = -°- 02382 


<p 2 l d| . i - i + 6 . 8. 6 1 J_ + 

2 7 5 13 160 8 250 


(3) 

f 

Jl/2 

+ 15 01158 

(4) 


| l/2 « - I - H • B - 150 + I V® - °- oklS3 
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(5) 


d (^. 2 I i 2 dl = /( - 6) (l 1 / 2 - 3S 2 + 35 7 / 2 - i 5 )i 2 as 

a 5 J 

- - + f i 5 - |fi 13/2 ♦ |t 8 


r 1 a(i£i e 2 al , - it + li . 36 + I . IB i . 3 _1 

Jl/2 d ^ 7 5 13 4 5 32 “ 


^fe + 33I5)-- 0 - 06683 


(6) 




(7) 


I 2 &* . 2 „ £ . — + -4/2 = -0.20013 
x /2 5 7 160 56' 


1 «) 1/2 -hi' 1 *’ '°- 68566 


Abbreviating 


U A ~ U 1 
u D - u-l 


to ( ) gives 


- b 1 ()()*A + b 1 »() 2 B - ^OO’C - 


+ b-,’() 2 F + — ^ bn'OG = K() 2 H 

1 ■ u o - ^ 


with A = 0.03374 

C = 0.02317 
F =- O.O6683 


B = - 0.02382 
D = 0.04783 
G =- 0.20013 


H = - O.68566 
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Introducing (45) and (46) 

tl/r ° = ; - 


172 


u o 1 


d,( ) + ^ 


*?■>] 

d(bi/r 0 ) = _ ( y 1 *l( ^ + | e l} 

d ( x / r o ) ( ) 3/2 [] 3 / 2 


gives 


( ) { } ( )2 b _ 1 


( ) l/2 [] l/2( } ( 3 A ( ) 3/2 [ ] 3/2 ' ' ( ) 1 / 2 [ l 1 / 2 


( )( )*C 


Un 


( )' { 1 , ”1 ( )' { } 


1 1 1 l ^ V V ) L J f “X V ) l J / \ r 

u o - “1 TW TW B 375 * ’ “o - u i Bp 71 


= K ( ) 2 H 


or 


( ) 

( ) l/2 [] 3/2 


A( ) 


uo - ui. 


B 

u o 


+ c( ) > ♦ £*] + 

^0 o J u o 


>1 . ui 

1 C 1 

+ B( )• 

n_Ji 1 

S' 

0 

I 

£ 

1 

0 

3 

Uo 


)♦ 


U r 


H i eX 

"o 2 "I/ 


^1 

Uo 


* ^ 1 4 * 1 -} 


= - K( )% 
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or 


( )’ 


r 


( ) 1 / 2 n 3 / 2 


- U 1 


*l( 


A + B + C + F 


>< > a ♦ £<*(• 


A + S + C + £]+ 


d i (D + G) > ( } + r < e i( D + §)> r' K( )% 


■A U.-1 u. 

o 1 o 


hence 


( ( )^ 2 n 3/2 1 


a 


«o " u l, ,2 


5 + + ^ 


Ul U-I 


"o - U 1 "o 


r = -K( )% 


with a =-0.004505, p =-0.01736, 7 =-0.01343, H =-0.68566 


and finally the differential equation 



= K d ( x / r °) (47) 

with oq_ = 0.00657, 3i = 0.02531, 7 1 = 0.01959, dj = O.13352, 
e ± = 0.25714 


Writing y instead of 
the equation 

- dy 


U A - U 1 

> ~ «i. 


and 


x in place of x/r Q involves 


Ay 2 + By + C 


y2(Ey + F) v^T 


Fy 


? = K dx 


or 


- dy 


^ Ay 2 + By + C 1 


Ey3 + Fy 2 /I J y 2 + F y 

» E > 


H K dx 
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vhere 


A = a]_ 


~ U 1 


u_ 


B = Pi 


Sk 


u n 

c = n — i 

"*■ 11 


u i 


Uq u q - u x 


E = d-, 


u o ~ U 1 
“o 


F = e x — 


^1 

"o 


Transformation of 


y = 


_ - f/e 


i - V 


= fl/y)/y 2 +'F/Ey 


dy = 2 F . t 


di E (i - £2 y 


y 2 = 




(i - 1 2 )- 


3 - (F/g)3 

(1 - t 2 ) 3 


gives 


Kx 


(i - i 2 ) (i - e 2 ) 
^ _ E (F/E)3 + F (f/e) 


2(F/E) 


(F/E)- 


(l - |2) 3 (l - l*f (l - f) 


(l - l 2 )‘ 


d| 


2 

a(f/e) 2 - b(f/e)i 

[i-i3 

1 + c(l 

- 1 2 ) 

Pj 

1 - e(f/e) 3 + f(f/e) 2 (i - I 2 ) 

i 


d| 


h\ 

2 

fS, 


f(f/e) : 


.(i _ I 2 ) + F ( F /E) : 


e(f/e) 3 - b(f/e) 


f(f/e) : 


A + 


(“ - B)l 


dl - 


- E(F/E) + F^l - £ 2 ) 


dl 



/CE _ 
\ F 


B 


f(f/e) 


dl - 
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A + 


*■ 2 1 

fE F 

(r 


2C - B(F/E) 


(f/e) ~ 


I + 


3 (f/e) 2 J 


? 3 + 


(f/e) . 


1/5 


r + const. 


By retaining the old variables and inserting the original values for 
the constant s, we obtain 


X / r ° K V uc 


u o - u l\ 1 ^ 2u o 




(tI/2 


£(D 


372 / 


+ x K /i 


( 48 ) 


with 


5 = 


U A - U 1 

k Uo - U1 


Ua - u n 


U D - u-L 


+ n 


U 1 


Uq - Ui 


5 ( 1 ) = 


1 + n 


U 1 


u o - U 1 


and the constants 
a = 0.0275 3 = 0.054-9 


7 = 0.0375 


H = 1.9259 


, U A " U 1 

The integration constant was so determined that the value I = 1 


. , "o " u l/ 

was obtained for the core boundary x K/ r o> x£/r 0 to be taken from the 
theory of the central zone. 


The formula still contains the empirical constant K which must 
be defined by experiments (section I). 

For the case of quiescent outside air ( — = 0 ) the differential 

) 

equation reduces to 


- di 


U A - U]\ J Oq 

u o - U l/L 3, 


d 1 3/2 ( S A " V 


\ u o - U 1 


= K d 


( 7 r °) 


■ 11 
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which gives 


s / r o -&A) 


F 2<< 


U A - V 
^ ~ u l; 


( x K/ r o) 


(49) 


with F 2 = O. 13 V 7 . 

(b) Temperature Field 

I. Equation of conservation of heat (44b) 


N’-of ^ + ^ 1 a 1 


= 1 ( 50 a) 


By (9) 


<pU) = 


t = (l - n 3y/2 ) 
(l-[o ^] 3 / 2 


2 0 


1 

/E 


\fE " 


( 50 b) 


Evaluation of integral: 

( 1 ) 


Jw dT) = J (1 - [. /k ] 3/2 ) 2 (1 - f i, *1 

=. £ . 5*3/4 7/2 + 1*3/2 5 . ^7/2 + 4*3/4 5 . 
27 5 7 ‘ 5 

isVV 3 / 2 + in 5 - ^3/\i3/2 + ^ 3 / 2,8 

/: ^ 43 = ^ 43 = 5 ( 9 / 7 °) + ^ (##«) 


E' 


5/2 V65 X 8 


III II I I ■■■■■!■■ Ill III 


111 ilium 11 1 mu 1 i iiihbih 
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For E = 2: 


Cp\|lT) dTi = 0.04382 


dtj = J ( 1 “ T i 3 ^ 2 ) 2 t dT i =. | n 2 - - n 7 / 2 + 5 t )5 


+ti dn = 9/70 = 0.12857* 


hence 


( b 2/ r o) 2 ( l3 A^o 

t£ / r ° ’ TFlTu 


u o “ U 1 / U A " U 1 


u o V u o “ u l, 


U 1 

dp + — ei 
Uq 


a. /a 

A/ o 


Uo - ui/% - ui\ ui 

u o V u o " u 3 2 u o 6l 


with d 2 = O.O 8765 


In addition 


e ± = 0.25714 


d ( b 2/ r o, 

d(x/r 0 ) 


= - 1/2 X 


"0 ~ U 1 , A [ U A “ U 1 \ , UA - u x ^- /jt ^ . Ul - b 

-sr~ wz^l l gzaiiLiZ -i ^ ^ A ' 

“o - U 1 , - u l\ U 1 1 3/2 I „ \ 3/2 

hr - ^ ei J W 


This value was incorrectly given as 0.25714 in the originial 
German version of this paper. 
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II. Heat equation (44a) 



We specialize to tj = l/2. 

Evaluation of integral: 

( 1 ) 

an £ (4 - [n a„ = - ^3/\ 7/2 + 1^5 

£ /2 <Pn a, - i - ■+ !^ 3/2 


For 


+ = 


E = 2: nl/2 1 x 

Jo tpi *> ’ 5 - 14 V s + to / 5 = °*°5773 

(i - n 3 /*) 2 - 1 - ^ - ,3 ; Va . | -j£. 


0. 41789 



= 0-02413 
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( 2 ) 


Jf^dr, = J( - 3 )(i - [ti/e] 3 /^ y/F^rj 

= -3E 3 ^J(rp / 2 - E 3 /\ 8 / 2 )dTi = - |e 3 /\ T / 2 + *k 3 / 2 Tp 

f l/2 JVdr, = 4/2 E 3 A + -1^3/2 

Jo d V 56 r 160 


For E = 2: 


/ ^’> 2a, > - • is V* + ^ * •°- 07438 

/ l/2 ^)*i/a- 0 - 03708 


(3) 

1 1/2 
= 0.04382 - 
Pl/2 


.//F 

eptT] dT] = / Cp4fT] dT] 


I',*" 1 ’- /i 

r 


I 


1/1/E Pl/2 

(pi 7 tj dt) - / cpiJ/Tj d^ 

Jo 


CpiT] dT) 


(for E = 2) 


11 4,1 ' (g - lik + ife) + ^ (&j - £3 + I ik) 
£5 - 3/ ‘ °‘ 04157 


r 

J 1/2 


cp\[rr] dr) = 0.00225 
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( 6 ) 


f -^E^- T^dT) =-0.023.67 

J 1/2 dT l 

3^/2 + 3^^^ = - | tj 7/ 2 + | n 5 
■^r + 1 " lis ) + J^y^=- o-2ooi3 



(7) 


Trp = - 3 ti 3 / 2 + 3tj 3 

dt] 


[&) -|- 1^.-0.68566 

^ d Vl/2 ^ 


Ua - Ui ' 

Abbreviating ( - — — — ) to ( ) gives then 


u “ U T 
. o 1 


- to fA ( )*A + bp’ fA ( ) B - 


( )f ^ + ( >(v°)' 

u o ' > * _ 


C - 


3 2^-- - u -(iA) ,p + b 2 ,( } ^ F + b 2*- — - ^ 0 =vi K( ) !a H 

“ u l\ ' > d n u o “ U 1 *o So 


with A = 0.02413 
D = 0.04783 


B =-0.03108 
F = - 0.02167 
H =- O.68566 


C = 0.00225 
=-0.20013 


or 


b 2 ( )' a + b 2* |A .( )P + b 2 (yo)'( )? + b 2 ~% - Ul /^oj 5 + 


b 2 ’-i*— — € = ( ) 

^ Uq - uj. d Q d 0 
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with a = 0,02638 3 = 0.05276 7 - 0.00225 

5 = O.OV 783 e = 0.20013 £ = O .96967 

The introduction of (51) and (52) 
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with a x = 0.00657 Px = 0.02531 = 0.01959 

= 0.13352 e 1 = 0.25714 d 2 = 0.08765 

Pl = 0.02413 P 2 = 0.05224 P 3 =-0.00199 

?x = 0.96967 



( 55 ) 
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1*0 


with the functions 



The integration constant was so defined that 

V« sl 

Hence, for infinitely great nozzle distances in the case of outside air 
in motion f u^ ^ 0 \ 
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b 2/ r o ~ \/? c 


By this result the assumption ^> 2/^1 = based upon the theory 
of asymptotic distribution functions, is sustained. 

But this does not hold true for the case of quiescent outside air 
( u i = °> The theory of asymptotic distribution functions gives an 

impracticable result: b 2/ b l~ * °°* The assumption of the same ratio for 

t> 2 / b l as for outside air in motion proves inconclusive. 

In the case of quiescent outside air, the differential equation 
reduces to 


, (V h ^ 

1 V A / 0 / 

“ u ! ' 


iU 0 - u i y 




U A - U 1 
u o - U 1 


with 


<Xi dp 

72 = 5l TW 


The solution reads 


\!*o 


U A ~ U 1 
, u o ~ u l, 


+ (* - 


For the breadth of the mixing region by 


b 2/r 0 = -= 


172 /u a - u^l/2 .1/2 


(w 2 
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and 



In the asymptote 







and 



With the constants computed on the assumption bg^b^ = \J2 the result . 
would he h 2 ^hj = 1.98* hence contradictory. The assumption ^ 2/^1 = \^2* 
for the case of outside air at rest must therefore he abandoned and the 
correct value of b 2 /b^ obtained by a special consideration. 


For this purpose the equation 


di 


1/2 


b2 / r ° = 7175 72 V° 


dg 


and ' 


b 2 /bi 


$1°! 

d lPl 


( 58 ) 


must be solved. The constants at the right-hand side contain b 2 |b 1 as 
a parameter. By an approximation method b 2 /b-j_ ~ 1-33 . 
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The calculation made with this value of gives for the case 

of outside air at rest: 

/u a - u l' 


00o - Ul 


v ‘‘ " W T 7(PP P 


b 2/ r o = 


with 7 = 3.245 


hence in the asymptote 


In addition 




/a a - V 

\ i/2 

7 1 

r 2 + (! - 7 2 ) 

(4 - 4 

] 


/% ‘ u i\ 




\ u o - u l/ 




7 2 = 1.123 


lj r o ~ !• 


33 b, r 


1 x o 


(59) 


(60) 


(61) 


a/\ 


0.79 


U A 

-M 

u 

- U-, / 

o 

1/ 


( 62 ) 


as against 


5 a/ s o ~ 


"A " U 1 


obtained for outside air in motion. 


2 \«o “ u l y 

The calculated functions are represented in figures 1 to 6. 

The integral in formula (55) was obtained by numerical integration. 


III. Calculation of the Asymptotic Distribution Functions 

At very great distances from the nozzle, the velocity and temperature 
differences of the jet and of the surrounding medium can be regarded as 
small. 

Limited to small temperature differences (density approximately 
constant) the differential equations (part 1, reference 1, (24)) read: 
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momentum: 


heat: 


mass: 


where 


_8u -8u / \ j 

u + v— = c(x) i 
ox or 


S 2 u 1 8u i 

dr 2 r Sr 


-bd , -bd - / \ j 
+ v— = Ee(x) > 


S 2 -a l 

Sr 2 r Sr 


c)u — Sv 
r — + v + r — = 0 

Sx Sr 


(x) = K b 1 (x)^u A (x) - u-l) 


(63a) 


(63b) 


(63c) 


1. The velocity distributions are now to be calculated, 
a. Outside air in motion u^ ^ 0 

The similarity formula 

ft - u i u A (x) - u ± 

^o “ U 1 “ ^o * U 1 

r 

with T] = 

v = v(x) x(n) 


b x (x) 



( 6 %) 


(64 b) 


gives 


d( u / u l) 

br 


A (x) - u x dcp 


U 1 
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^(u/u-jJ u Q - u-l u A (x) - u x d 2 cp 1 

— — ^ 

or uj Uq - u^ dr; b]_ 

c)(u/ui) _ u Q - u x / 5 a (x) - ^ , u Q - / u A (x) - u-A dcp / V \ 

8x U! \ Uq - ui / ui \ Uq - u x /dT) \ !)]_ / 

etc. 

On the basis of computing the axial functions by the momentum and 
the heat equation 


U A (x) - 2/3 


~ x 


b-, (x) ~ x 


1/3 


u o _ U 1 • 

is applicable to the asymptote according to (15) and (l6), 
By the equation of continuity of mass 

a A 


H fi - u-jY / u a u lY ' .. v(x) 


"o " u l> 


, u o - u a 


Ui 


hence 


iki ~ X -V 3 

u l 


for the asymptote. 

By the equation of motion 


U A 

- Ul \ 

u o 



/S_ 

- 

I*!* 

l u o 

-u 1 ) 

Ibl 


« -/ 


y— 



/— 

\ 

/ U A 

- u-^ 

1, /. I 

< U A 

- u l\ 

\ u o 

- u iy 

)i/ b i 

V u o 

-uj 


( 65 ) 


The first two terms and the last term are of the order of magni- 
tude x-5/3, while the term —|l/bi is of the order of x - ^/^ 

U 1 \k> ~ u l/ 

The term with v(x) is therefore disregarded in asymptotic con- 
siderations. 
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Noting further that 


u - Ui \ u Q - u-L 


*1 V u o - u l/ U 1 


+ 1 a? 1 


for small velocity differences gives for the case of outside air in 
motion the asymptotic motion equation 

- "i kfcfa) +fl/r) feil. 


the relation 


u A(x) - u x _ 2/3 


u o - U 1 


h x (x) = h M x 2 /3 


gives 


, 1 -2/3 

2 -5/3 -2/3 Do& 3 x ^9 

- f u m x J q) - u x ' ~j. 5 t] ~- 

3 00 “ b xV3 drj 

oo 


T , u o - U 1 a 2 -V3 ] d 9 , 1 d <P 

— : — , ; U . X < — > — 


1 booX 


r, 2 I dT l 


9 1 d.cp Uo " U 1 1 „ J d 2 tp 1 dcp 

_£ C p_ r , K — U 00“\ — p + 

3 3 dr, u-L h OT 1 diq2 ^ dr| 




\ ^Uo - up 

r- K ^r--~r 

— 


i| + ^ for r, ^ 0 

dr, 2 dr. 


Integrated 


1 2„ TA Uo ~ U 1 

r) qp = K 

3 Ur 


— n — 

h m dr 


i -a. -q dr, = 2s 

3 u „ K u o - U 1 9 
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Hence 

cp = e**( aoTi ) 
But by (15) and (l6) 


2 t t u-i b 

with cr_ = — — i Si 

0 ... 6 K u o " U 1 


1/3 


(67) 


te) 

(Y^f 75 


- kV 3 \ y J oh with ^ = k.2h6 


ji/3 

i?73 /u 0 - u-^ 2 / 3 1 


1 ( u l/ u o)‘ 

U co = ^ " o ' /o ' E 2 Vlth B 2 = 0.216 


u r 


Therefore 


a Q = 1.811 


(68) 


Transforming correspondingly b-^ 
gives 


xV 3 


on the variable T}* 


= r/xl/3 


9 


= e’^ 0 *^ with cr 0 * = 0.427-^- . (V U °) 

r1/3 ^ ~ Ui j!/3 


(69) 


b. Outside air at rest (uj = o) 
Conformably to (28) and (29) 


u A ~ l/x 


b^ ~ x 


we put 


u = u A (x)q>(T)) v = v( x)x( T] ) with T) = r/x 


(70) 
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A similar consideration as given under (a) yields 


v(x) ~ l/x 


(71) 


The apparent kinematic viscosity follows as 


e(x) = Kb-^(x)u^(x) = constant 


Thus the problem appears to be reduced to that of computing the 
distribution function for the laminar viscosity e. This has been 
solved by Schlichting (reference 9) • 

The solution is 




But according to 

b ro = K 20.321 

hence 

a Q = i 0.0784 (73) 



2. Relative to the temperature distributions Reichardt (reference 6) 
has shown that 

\|r = cp Ar / A Q (74) 

with A T , A c being exchange quantities of momentum and of heat, 
respectively. 

To put it briefly: In this representation, the ratio A^JA t 

corresponds to the factor E. With E = 2 

* = 9 1 / 2 (75) 

The calculated distribution functions are represented in figures 7 
and 8. 
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B: GREATER TEMPERATURE DIFFERENCES BETWEEN <IEr 

AND SURROUNDING MEDIUM 
I. Methods and Results 


I n view of part A of the present investigation, which dealt with 
the diffusion of a hot air jet for small temperature differences, we 
can he brief in many respects. 


The differential equations of the mixing field read: 
equation (2)) 

equation of continuity of mass 

S(rW + a( rev). = Be(x) |_J r 3p. 

ctoc 8r ° r 8r 

V. J 

equation of continuity of momentum 


8(rpuu) + 8(rpvu) = 


8x 8r 

equation of continuity of heat 

8(r£>uT) S(rpvT) 


dr 


_du __ 
rp — + Eru- 

dr drj 




dx dr 

with the apparent kinematic viscosity 


= Ee(x) — It- 


dr 


_d(p!F) 


dr 


(part A, 


e(x) = Kb-^x) 


u 


max 



Integration with respect to r gives for the transition zone (part A, 
equation (3)) the moment-urn equation 


rpv 


( a - u i) - 55 

the heat equation 


a f b i „ 


pu^u - u-j^r dr = e(x)r 


d(u - u n ) 


dr 


rpvd - — 
dx 


32 8 

pu£r dr = Ee(x)r — ( p3 ) 

or 
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with 


e(x) = Kb 1 (u A - 


For r = 0 a further integration (part A, equation (4)) gives the 
equation of the conservation of momentum 



dr = 




and of the heat 



pu£r dr = PoUq-Sq 



The investigation of the mixing field was then carried out for 
small temperature differences (part A) by means of the cited integral 
equations. For the velocity and temperature distribution the "similarity" 
formula was used: 


* ~ u i = ~ u i y (T)) 

u o- U 1 u o - U 1 ^ n 




n = 


r 



This formula expresses the following facts: For very great nozzle 

distances (asymptotic case) the flow proves itself similar according to 
theory and experience. But it also proves itself almost similar for 
nozzle distances extending up to a point near the boundary of the core, 
according to experimental data (cf. Pabst, reference 2), 

In order to simplify the calculation, the asymptotic ratio 


b^/bi = v/eT 


for the total transition zone was assumed. 
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The functions m and i|r 
tribution functions (part A) 


are identical with the asymptotic dis- 




n 



The application of the foregoing investigation method to fields 
with greater temperature differences involves a considerable amount of 
paperwork. For this reason, a different method is applied. 

The study is limited to fields with outside velocity u^ ^ 0 
different from zero, while the singular case u^ = 0 (cf. part A) is 
disregarded. 


The goal is to calculate the breadth of the mixing regions of 
velocity bj and of the temperature bg as well as the axial 


functions 


u A (x) - Uj 

U — U-. 

o 1 


and 


*A( X ) 




For the calculation, the differ- 


ential equations are used which for the jet center read 
_ dp(u - ux) ^ 1 


u- 


8x 


€(x)2 


d 2 p 


o (u - Ui) 

P-1 + E(u - u x ) — 


8r‘ 


8r c 


u|^ = Ec(x)2^ 


P 8? 




8r £ 


with 


e(x) = KbjXx) (u A (x) - u-jj 

also, the conservation formulas 
momentum 

1 r 'c 

/ pu(u -u 1 jrdr = P o u o(u 0 " u l ) ~f 


( 1 ) 


( 2 ) 


J q 2 ~PV$r dr = PqUq'Aq t q z Jz 


and heat 
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In conformity vith the almQst similar behavior of the flow in the 
transition zone described above the following proposition is made: 


u - u-l = (u A (x) - u^ qp( q); t) = r/b-j^ 

(3) 

d = ^ A (x)qp(ii*); ry* = b 2 


In other words, it is assumed that in the whole transition zone 
the distributions over the breadth of the mixing regions are represented 
by the asymptotic distribution function 


rC^oTi) 2 and e “ (^o^*y 




The constant cr D can be determined by experimental calibration. ^ 
Approximating the experimental distributions by the function (l - t] 3/2) 
gives (cf. part A). 


a = 1.8l 

With (3) the following differential equations 

(\ -”iY 


Uo - U X 



and 


a ( x / r o) 

(*A - u l' 


lu„ - u, 

\ o 1/ 


1 

K 


1 

2 



U A ~ V 
u o - 


U-, 


u o - "l 


U A - U 1 

Uo - Ui 


are obtained from (l). 


d^cp 

dT] 2 


Tj=0 


(5) 


( 6 ) 


The first equation defines the mutual correlation of the axial 

functions |— — ] and second describes the variation 

f *A - u l\ 


/U A - u l\ 
v u o " u l; 


of the axial function 


^ - U l; 


with the nozzle distance x jv 
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The conservation formulas (2) give by reason of (3) 



with the coefficients 



with the coefficients 



(7b) 

[i] = 2 r 2di dn* 

Jo [1 ^o/T^A/"o)cp] 

^ where , with tj* = r /b 2 

t(T]*) = cp^rv* b 2 ^b 1 ^ = qp(rj)^ 

In complete generality, the solution of the described equation 
system (6) and (7) still presents a difficult problem, which is largely 
due to the calculation of the integrals by which the coefficients of 
the breadth of the mixing regions are represented. 
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By contrast, the solution for the asymptote (very great nozzle 
distances) is readily indicated. 


With 



The differential equations (6) give then immediately 



With the nozzle distance x/r 0 as independent variable 
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(10 con*t) 


Consequently, greater temperature differences for the asymptote 
make themselves felt "by the factor rTs* It implies that 

(i + 

the breadth of the mixing regions increases slower at greater tempera- 
ture differences and that velocity and temperature decrease considerably 
more. 


The asymptotic solution is at the same time to be regarded as first 
approximation for the general solution. Comparisons with better approx- 
imations indicate, however, — for example, in the theory of small temper- 
ature differences (cf. part A) — that the asymptotic solution still 
represents no satisfactory approximation. 

To gain a second approximation, it is logical to continue with the 
equation system (6) and (7) on the approximate assumption that 

b 2 /bi = \/E (which holds rigorously for asymptotic conditions). 


For E = 2 — this value is obtained by experiments as will be 
shown later — the coefficients of the conservation formulas can be 
indicated analytically: . 


(id 

(i) 

•M 

H 


' 

2 < _ 

3 (V^ 

J 3 - i (v t i: 

) 2 + ( 5 a/ T i) - W 1 + 5 aAi) . 

OJ Al 

O 

to 

' 5 a/ t i) ‘ 

(®- 

- Jn(l + I A / T ij 

J 

2n(l + \/ T l) 

“° 2 1 
1 

( ; 

i N T h 

V T i) 2 

i 2 - (\hi) * 

°o 2 

1 l ± 
°o 2 

1 ( 1 + ^Al 

(Vi 

(Vif 

Til 

) 


Introduction of the functions thus obtained for b-^ and b 2 in (6) 

gives a system of ordinary differential equations which, however, can be 
solved only approximately. 
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But a second approximation in analytical form can be obtained, if 
one limits oneself to the asymptotic coefficients of the conservation 
formulas which read: 


(II) =2 <p 2 n dri = \ \ 

JO 2 cJq 

n°° 

(D = 2 / cprj dri = l/a Q 2 
n°° 

|"lll = 2 / q^T]* dt)* = i 

1 A JO 3 o 0 2 


( 11 ) 


[i] = 2 J Q w* = i/° 0 2 

This solution seems to represent an adequate approximation. We get 


h ll r o 


(1 + *JT 2 ) 1/ ' 2 f u A - u l \ l/: 
1 V^o " U 1 


■- (ii) + — (i) 
u o \ u o " U 1 / u o 


'll/ 2 


with (II) = - 1 


2 a 2' 


(I) = l/a c 


^ i To = 1 1 1 

2 I ° (l + * 0 / T i ) 1 / 2 (^a/^o ) 1 / 2 r^o - U x /n A - Ul \ rn + Ui 


( 12 ) 


Ur, \Un - U 




1 / 


172 


“ ith N - | ^ 


[1]- 


With these functions, the differential equations (6) give: For 


the correlation of 75 . /t3 N 

\ A I°; 

(*.l* \ - 1 

) and " U l' 

) 

L 

\ A l °) - E . , 

/ _ 1 - l\ 

, + 1 - "hn 

1 ^o - u iy 

if 1 + 2 Uo - Ul ) 


( U A . U 1 
V \ u o ~ u l, 


6 


/ U A ~ U 1 
\ u o - ^l 


+ 2 


U 1 


u Q - u x 




( 13 ) 
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ip 


(V °)* 


The integration constant was therefore so determined that 

/ U A ~ V 

\ U 0 “ U l; 


= 1 for 


= 1. 


The dependence of the axial function 
distance xjr Q reads 


U A - U 1 

Uq - ui 


on the nozzle 


7 r o = 


[l + *opl) 


1/2 


1 71 jl ~ 

K 8 a 0 V Uq 


- -l\ l/2 ‘ 


) -fe -5(D 1/2 ) 



(i4) 


where 


5 = 



for 


The integration constant being so 



(boundary of core) 


defined that 


x/r Q = x 



The two empirical constants K and E appear in the theory; 

K = 0.010 for asymptotic conditions (cf. part A). The constant E 
equals 2 according to the experiment. 


The structure of the mixing field for 


0.75; 0.95 
0.75; 1.5 

obtained with these constants is represented in figures 13 to 25. 

On experimental data the measurements by Pabst (reference 2) are 
available. He measured the diffusion of a hot air jet of s; 300° C in 
the transition zone and a relative speed of ~ 4oo m/sec for l8, 101, 
and 107 m/sec outside velocity - 

l.O; "° ~ Ul ^ 0.95; 0.74; 0.53 

I ^ 


"o - U 1 


= 0.5; 


Uo 

V T i = 0; 


■■■■imi 1 1 in i ■■■! uni i in iiHiBi iii i iii iii 
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Pabst’s measurements of the velocity and temperature distributions 
over the breadth of the mixing regions were very complete. His measure- 
ments at the last three sections (x/D = 16 ,. 20, and 24) are reproduced 
in figures 26 to 28, to which we have added the theoretical distribution 
functions for the asymptote 


and 


qp = e T) ~ r/bj 

* = qpl/E (E = 2); 


a being defined accordingly. It is seen that the theoretical distribu- 
tion functions reproduce the experimental distributions quite closely. 


The significant result 


E = 2 


is quite plain. 

involves mainly a compar- 
f unctions 


Pabst ' s average values are represented in figure 29. Figures 30 to 
32 contain the comparison of the theoretical and experimental decrease 
in velocity along the jet axiB, the constant x - K/r 0 of the theoret- 
ical curves being defined accordingly. Surprisingly the measurements on 
the whole indicate a greater decrease in velocity, although this was to 
be expected since, as a result of the friction, at the inner and outer 
nozzle wall some loss of flow must be reckoned with, which stipulates ap 
effective radius different from the geometric radius. Figures 33 to 35 
show the comparison of the theoretical and experimental temperature drop 
along the jet axis, with the constant of the theoretical curve again 
properly defined. The reservations regarding any accidental loss on 
heat flow disappear. Even so the measurements indicate a fundamental 
departure in the sense that the experimental decrease is flatter; however 
this difference should raise no objection to the theory since the tempera 
ture measurements along the jet axis seem unfortunately to be faulty. 

This is seen fairly plainly in figupe 29 where, with increasing distance 
from the nozzle, the curves of the velocity and temperature drop approach 
one another and then even intersect. By contrast, the experimental dis- 
tributions give E = 2 somewhat plainly, which, in other words, means 


The check of the theory on the experiment 
ison of the theoretical and the measured axial 

u.U) - u-, _ 

— and 

Uq - ui / 
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that the temperature exchange is greater than the velocity exchange 


(according to theory ^ ^a/^oJ 
nozzle distances). 


~ 1 r*A ~ U 1 

"o " U 1 


£ 


is to he expected for greater 


To he sure the described behavior is due, to some extent, to the 
friction losses hut not enough, according to preliminary calculations, 
to explain this difference. 


II. Calculation of Mixing Region Structure 

1. The differential equations describing the diffusion of a hot air 
Jet in the transition zone read (cf. part A): 


drpu(u - u x ) + drpv(u - u-J = ^ d_ 


dx 


dr 


dr 

drpufl drpv-fl 
dx dr 

dru + drv _ Q 
dx dr 




= Ee(x) |p -j r 


dr 
d(pB) 


dr 


vith 


e(x) = Kb^x) (u A (x) - u^ 


or 

ru 


_ dp(u - ux) _ dpfu - up) , . d _ d(u - ux) , x dp I 

u — i ±L + rv — i L - e (x) J rp — * =*- + r( u - u-, ) — > 

dr dr ' ■'dr 


dx 


dr 


ru 


dx 


dp(u - u x ) 


dx 


dx 


+ rv =.Ee(a 

:) 

r ^ 

dr 

dr 

dr J 

= 0!) 



lim i £L_ < 

_ ^(u “ U 1 ) 
rp — - 

r _> o r dr 


dr 

= E€(x) lim 

l 2_ . 

dp3 

r 

r — > 0 

r dr 

dr 


dp 

'd7 
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Considering that the first derivatives with respect to r disappear 
for the jet center, one obtains 


_ 8p(u - ux) = 6 
8x 


(x)2 .p!JlpJ + E(u-u 1 )!-§. 


with 


If one puts 


_ 8p$ -r, / \r, - 8^ - 

U = E e(x)2-S p yy + i3 rr 

a x 8? -a? 


e(x) = Kb 1 (x) (u A (x) - 


— _ const. 1 

P " T x 1 + fl/T-L 


according to p X T = constant, the equations read 


u - ux 

d l * 3/Tj . . 1 1 ag ( 5 - 

ax W 1 ♦ a/T-L ar2 


!(u - u-l)- 


.2 1 


U = Ee(x)2-j 

8x 


5 2f 2 

VI + */?!, 


For computing the mixing field, the conservation formulas are 
employed again: 


momentum 


pup - u x )r dr = P 0 u 0 (u 0 - u-^ r 0 ^2 


pudr dr = P 0 u 0 ^ 0 v 2 j 2 
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1 - ”l) r ar , I "lJ ro 2/ 2 

Jo 1 + V T 1 1 + *o/ T l 


U-a « 0*0 p/ 

— — -r dr = ■r Q e y 2 

to i + , ®/t 1 1 + ^q/T-l ' 


2. The lav of similarity is applied 


“ U 1 = ( U A^ X V - u l)<PU)> t) = r / /h 1 (x) 
■5 = ^ A (x)cp( ry *) , T|* = r/b 2 (x) 

* = \(x)q^T] = \(x)t(n) 


hereby it is to be noted that 


(8/8x) r = (s/ Sx )t) " (V a i))n~ 

(8/8rJ x = (S/8T))i- 

For the differential equations (l6) ve get 


<(u A -»,)»♦ Ul ) > +&/*» 


u A - UjJcp 


( U A - u l) 2 j 1 + (^ T] y 


' t A 

8^/u a - u,)<p , h +( a Jf-Ay 


= EKb 


l( u a - u l) : 


f *( y. .a 
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or, bearing in mind that ,<p = i|r = 1 at the jet center 
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dijr _ d<p(ri ^ 1 /^ 2 ) _ 1b^ dtp 

^2 d^ 


dtj 


dT^ 

d 2 cp _ / "bl \ 2 d 2 cp 


d n £ 


, b 2/ drf 


we get 


u n 




( 

' U A 
Uo 

- U. ± \ 

- u iy 

t 

1 ( 

‘*A ■ 
^ • 

t 1 

</CV°) 

[ 1 + 


] 

[1.1 

[y^ 

(V 

M 

-2 


V. 


d 2 cp 

dt)2 


n=0 


E 


f^A ~ u i\ (•' 

>oAJ(V*° 

) 

v^o - u v 



M 

r 


1 d 2 cp 

2 


T )=0 


corresponding to 


/, / u A - u l \ “1 


\V 


= 2EKb 1 ( ' UA 

1 \ u o ~ U 1 

(VT])(?a/^) 


\ 


( ^a/^o) ’ 

• -i/[[i +( 

( ^ a/^o 

,v T j)(fy 4 °2 

) IS 

L 1 + ( 

V Tl )( 

[A/ A)] d 


K/ (, °)( < W T i)( ? a/' 8 c 

>)’i 

[X + (- 

3o/Tlj(-T A /^o) _ 

2 

_ 


T )=0 


1 d : 


dm 


[1 + (V T i(^A/%)] 2 b 2 2 dT) 


T )=0 


(The accents (’) signify derivatives with respect to x. ) 
Hence 



(20a) 
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- cw 1 ft>l\ 2 / *A ~ U 1 

b 2\W V* ~ U L 



T)=0 


(20b) 


Dividing the first equation by the second gives 


(|jj[ 

(®A/ fl o) [l +(S 0 / T ^a/ i3 o)_ 



finally 





( 21 ) 


The elimination of ^a/^o) from equation (19b) by means of this equation 
leaves 




( 22 ) 


Equation (2l) gives the mutual correlation of 


(^A j^oj 


U A 

- «i\ 

u o 

- 


and 


along the jet axis. 

Equation (22) contains the dependence of the axial function 
on the nozzle distance. 


U A 

- M 

u o 

- u i/ 
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3. The application of the equation of similarity (19) to the 
equatiops of conservation (l8) gives: 


momentum 


similitude 


1 a fe- i ;i). rar , M u °- . ' i i) r o 2 / 2 

0 1 + *o/Ti 1 + 1?o yT 1 1 


u - Ul = (u A - uj cp(tj), t] = r/h-L 
= a.(x)i|r(T)) 


It yields 


r = ^(Up- u xi , 

Jo 1 +( i3 o/T])( 3 a/ i3 o)^ q ^ 1 1 + « o yti / 


( b l/ r o) 


2 / U A - U 1 


^o ~ 


u i/l / [i +( < 0 /t^C« a /* 0 )^| 


dr] + 


u i / u a - u i 


^o - u l\ u o - u - 


&T)\ = 


1 1 


"l + ( l3 o/ T ^(\/ a o) ltr ] ) - Ul 1 + ^o/ T l 2 


finally 


A. , \2 f^A U 1 

W r °)k^ 


2l u A ~ u l\ u o ‘ U l/ U A " U 1 


"o V^o " U 1 


[! + (^/ T ^A/^o)^] 


dtl + 


ui / / \ = i i 

M / [i + (VTa)( V 6 o)>] / 2 1 + V T i 
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From it follows 
b 


Ir = 1 1 1 

11 ' (1 + V T i) 1/2 / °A - yi \i/2 rug - ^/gt - u x \ | ^ ,172 

v ' v^o - u i/ L u o \ u o - u iy u o 


with the coefficients 


(II) = 2 


(I) = 2 


9 2 n 


l 0 [1 


qpr) 


dii 


dfi 


1 0 [ 1+ '(V T ^V*o)*] 

The asymptotic coefficients are 


(II) = 2 / cp^T) dTj; (I) = 2 / qprj dTj 
JO JO 


with 


Hence 


heat 


<P 


. e -( c, ° r| ) : 


f 11 ) -1^5! (I) - ^ 


2 CTO 2 


°o 2 


1 + 


u-0 j ^^o 2/o 

__ r dr = - — : — r o y 2 


fl/Tl 


i + -a 


o/ 1 ] 


similitude 


with 


•3= ^(x)<p(t]*)j t)* = rj b 2 

U = Ui = ^u A (x) - u^tKti*) 

Sr ( tj*) = “ <PU) 


(23) 
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The result is 


r < ? . v° 2 /a 

/o [1 TW A >o) *] 2 " ' 1 - ^o/T! ° ' 


or 


£N*°)(f o [i + (v^VMk] 


dry* J = i = - 

/ 2 1 + ^o/Ti 


Hence 


b o/ r o = 


(l + 3 o/Ti) l/2 (\/%) l/2 


Uq - Up/u 

A - upN , 

ci] + —Til 

l/2 

Uq \U 

0 - U J ^ 

J UqL J 



with the coefficients 


fill = 2 f - dtf 

L J JO [l +(V T l)(^A/ l3 o)t] 


(24) 




94* 


/ 0 [1 +(V T ^V°)*] 

The asymptotic coefficients are 


diy* 


with 


which with E = 2 give 


jjtlj =2 J 1 cp^Ty* diy*; jjlJ = 2 J CpT}* diy* 


cp = e“( a o^) 

t = e “ E C a on *) 2 

&- k 3 & 
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Asymptotic solution (first approximation) 


The asymptote (very great nozzle distances) follows from (23) and 
(2k) as 

xV 2 


bl/r ° ~ (1 + VT^l/2 ui) V"o 


D 2 /r 0 ~ 


Hence 


1 1 /M 

(TTl^p (v^p W 

/“A - u l'\ l/2 

ta/bi . W - U V , 

' (<a/»o)V2 

Insertion in the differential equation (2l) gives 


1/2 


l/°c 


u fl - Up\2 


U A - V 


."o - u l/ = 1 \ Up - U 1 


( V °)’ * ( V * o ) : 


or 



1 

/“A 

- «l\ 

N 

” E 




From (22) follows further 

d ( x / r o) 1 1 K\ l/2 i 

d / d A - ui \ 2 K \u a/ a D 

W - u l / 


U 1 


(25) 


(26) 


^o 


- U n 


1 - d o/Tl) 17 ^ " , M 

' - U 1 j 


37^ d^ 

dq 2 


Tl=° 


or 


x / r o = 


111 
K 3 a, 


o d 2 cp 

dq2 


1 f — U1 -1 

(l + %/T;l) 1/ 2 W/ - Up / u A - Up 

n V U ° “ U l - 

T]=0 X 
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5. (Second approximation) 

With the asymptotic coefficients the result is 



With these functions the differential equation (21) reads 



I 



TO 
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Separation of the variables produces 



/ U A “ u l\ 

When the integration constant is so defined that = 1 for 

\ u o " U 1 j 

(*a/*o) = 1 then 



(29) 

The differential equation (22) gives 



Hi W 
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hence 



Transforming 
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- Temperature drop along jet axis plotted against the vel< 

along the jet axis. 





Figure 6.- Breadth of mixing regions of the velocity and the temperature 
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Reichardts measurement (reference 7) 
( Outside a/r at rest (u t =0)) 


Theory 


K- 0.0/0 


• x ID-- 133 
ox/D^/7.3 
a x/D-30 
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7j=r/x 

Figure 9.- Comparison of theoretical and experimental distribution function. 


• ftuderis fest cur ye (reference 8) 
(Outside a/r at rest fit, =0)) 


K= 0.0 iO 


Figure 10.- Comparison of theoretical and experimental velocity decrease 

along jet axis. 
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• Ruden's test curve (reference 8 ) 
(Outside air at rest (u, = 0 )) 1 


K= 0.0/0 


Figure 11.- Comparison of theoretical and experimental temperature drop 

along jet axis. 


Ruden’s test curve (reference 8) 
( Outside air at rest (u, =0J) 1 


Figure 12.- Comparison of theoretical and experimental relationship between 

temperature drop and velocity decrease. 
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Figure 15.- Velocity decrease along jet axis for greater density differences. 



Figure 16.- Velocity decrease along jet axis for greater density differences. 
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Figure 29.- Decrease of velocity and temperature along the jet axis. 
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Figure 30.- Decrease of velocity along the jet axis 
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Figure 31a.- Comparison of theoretical and experimental decrease in velocity 
along the jet axis for effective nozzle diameter. 


a Pabst’s measurement, ( Reference 2) 
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Figure 34.- Temperature drop along jet axis 
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